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ABSTRACT 

We present the spectral analysis of three Suzaku XIS observations of 3C 111 requested 
to monitor the predicted variability of its ultra-fast outflow on ~7 days time-scales. 
We detect an ionized iron emission line in the first observation and a blue-shifted ab- 
sorption line in the second, when the flux is ^30% higher. The location of the material 
is constrained at <0.006 pc from the variability. Detailed modelling support an iden- 
tification with ionized reflection off the accretion disc at ~20-100r g from the black 
hole and a highly ionized and massive ultra- fast outflow with velocity ~0.1c, respec- 
tively. The outflow is most probably accelerated by radiation pressure, but additional 
magnetic thrust can not be excluded. The measured high outflow rate and mechanical 
energy support the claims that disc outflows may have a significant feedback role. This 
work provides the first direct evidence for an accretion disc-outflow connection in a 
radio-loud AGN, possibly linked also to the jet activity. 
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1 INTRODUCTION 

There are several indirect pieces of observational evidence 
that outflows /jets are coupled to accretion discs in black 
hole accreting systems, from Galactic to extragalactic sizes. 
Recently, the most direct evidence came from the Galactic 
microquasar GRS 1915+105 (Neilsen & Lee 2009), where 
the appearance/disappearance of a relativistic disc line and 
blue-shifted Fe K absorption lines are connected to the state 
of the source and the radio jet. In radio-loud AGNs, rela- 
tivistic jets are routinely observed at radio, optical and X- 
rays. But the first discovery of disc outflows in this class 
of sources is due to Suzaku. A systematic analysis of the 
4-10 keV spectra of five bright Broad-Line Radio Galaxies 
(BLRGs) showed significant blue-shifted Fe XXV/XXVI K- 
shell absorption lines at E>7 keV in three sources (Tombesi 
et al. 2010a), namely 3C 111, 3C 120, 3C 390.3. They imply 
an origin from highly ionized and massive gas outflowing 
with mildly relativistic velocities of u ollt ~0.1c. 

These characteristics are very similar to those of the 
Ultra-Fast Outflows (UFOs) previously observed in Seyferts 
and quasars (e.g., Chartas et al. 2002, 2003; Pounds et 
al. 2003; Braito et al. 2007; Cappi et al. 2009; Reeves et 
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al. 2009). In particular, a uniform and systematic search 
for blue-shifted Fe K absorption lines in a large sample 
of Seyferts observed with XMM-Newton was performed by 
Tombesi et al. (2010b). This allowed the authors to assess 
their global significance and high detection fraction of >40%. 
Tombesi et al. (2011) then performed a photo-ionization 
modelling of these absorbers and derived the distribution of 
the main physical parameters. The outflow velocity is in the 
range ~0. 03-0. 3c, with a peak and mean value at ~0.14c, 
the ionization is in the range log£~3-6 erg s _1 cm, with 
a mean value of ~4.2 erg s _1 cm, and the column density 
is m the interval A w ~10 22 -10 24 cm , with a mean value 
of ~10 23 cm -2 . The mass outflow rate of these UFOs can 
be comparable to the accretion rate and their kinetic power 
can correspond to fractions of the bolometric luminosity and 
is comparable to the typical jet power for the radio-loud 
sources. Therefore, they may play a significant role in the 
AGN cosmological feedback (e.g., King 2010). 

Theoretically, the complex coupling between radiation, 
magnetic fields and matter that should be considered to 
properly explain the formation of outflows/ winds from ac- 
cretion discs has not been accurately solved yet. However, 
simulations show that disc outflows are ubiquitously pro- 
duced and can be accelerated to velocities >0.1c by radia- 
tion and/or magnetic forces (King & Fabian 2003; Proga & 
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Figure 1. Combined XIS 03 light curves of the three observations 
in the 4—10 keV band, binned to the satellite orbital period of 
5760s. 

Kallman 2004; Oshuga et al. 2009; Fukumura et al. 2010; 
Sim et al. 2010). The detection of UFOs in BLRGs repre- 
sents an important step for models of jet formation and for 
our understanding of the jet-disc connection. From a theo- 
retical perspective, disc outflows are a necessary, although 
not sufficient, ingredient for jet formation. For instance, in 
the magnetic tower jet simulations there is the possibility 
that some plasma is trapped and dragged upward by mag- 
netic field lines and winds/outflows may provide the external 
pressure needed to effectively collimate the jet (e.g., Kato et 
al. 2004; McKinney 2006). Thus, models attempting to ex- 
plain the link between the jet and the accretion process will 
have to take these components into account. 

The bright BLRG 3C 111 was proposed by us in the 
Suzaku G05 as the best candidate for a follow-up study of 
the predicted ~7 days time-scale variability of its UFO, de- 
tected by Tombesi et al. (2010a) in a previous 2008 observa- 
tion. 3C 111 is a powerful and X-ray bright (2-10 keV lumi- 
nosity ~l-4x 10 44 erg s" 1 ) FRII radio galaxy at 2=0.0485, 
with a radio-loudness log-R =2.35 (Sikora et al. 2007). Its 
AGN is powered by accretion onto a super-massive black 
hole (SMBH) with mas£] ~(2 - 30) x 10 s M Q (Chatterjee 
et al. 2011 and Marchesini et al. 2004, hereafter Cll and 
M04, respectively). It is important to note that 3C 111 has 
been recently the subject of an extensive monitoring cam- 
paign to study its accretion disc-jet connection (Chatterjee 
et al. 2011). Indeed, major X-ray dips in the light curve are 
followed by ejections of bright superluminal knots in the ra- 
dio jet. A similar campaign showed the same behaviour for 
3C 120 (Chatterjee et al. 2009). Both these sources show the 
presence of UFOs in the X-rays (Tombesi et al. 2010a; Ballo 
et al. 2011). 



2 DATA ANALYSIS AND RESULTS 

The details of the three Suzaku observations of 3C 111 are 
summarized in Table 1. They are spaced by ~7 days and 
we refer to them as Obsl, Obs2 and Obs3, respectively. The 
data were processed using version 2 of the Suzaku pipeline 

1 Chatterjee et al. (2011) used He? measurements. Marchesini et 
al. (2004) assumed the bulge luminosity relation. The discrepancy 
is probably manily due to the different extinction adopted. We 
suppose the former is probably more reliable. 



and the light curves and spectra were extracted following the 
standard procedursQ- Source spectra were extracted from 
circles of 3' radius centered on 3C 111, while the background 
ones were selected from circles of the same size offset from 
the source. In this letter we focus on the spectral analy- 
sis of the X-ray Imaging Spectrometer (XIS) in the Fe K 
band, E=3.5-10.5 keV. The analysis was carried out using 
the heasoft v. 6.5.1 package and XSPEC v. 12.6.0. The spec- 
tra were binned to a minimum of 25 counts to apply the 
X 2 minimization in the model fitting. Errors are at the la 
level for one interesting parameter and the line energies are 
at rest-frame, if not otherwise stated. The front illuminated 
XIS and XIS 3 spectra were combined after checking that 
the continuum and flux agree within 2% (hereafter XIS 03) . 
The data from the back illuminated XIS 1 chip were checked 
separately only for consistency. The combined XIS 03 light 
curves are reported in Fig. 1. We note a ~30% flux variabil- 
ity between Obsl and Obs2. 

A power-law continuum with r~1.65-1.7 and a nar- 
row (unresolved) Fe Ka emission line at 6.4 keV provide a 
good representation of the three spectra (e.g., Reynolds et 
al. 1998). The narrow emission line is most probably asso- 
ciated to reflection off material located in the outer accre- 
tion disc or the broad line region (Chatterjee et al. 2011; 
Ballo et al. 2011). We include a Galactic absorption fixed to 
N H =3 x 10 21 cm" 2 (Kalberla et al. 2005). We refer to these 
as the baseline models and the best-fitting parameters are 
reported in Table 1. From panels 1 and 2 of Fig. 2 we can 
note the presence of additional emission/absorption lines. 
This is confirmed in panel 3 of Fig. 2 by the contour plots 
derived with respect to the baseline models (e.g., Tombesi et 
al. 2010a). In particular, the contours suggest the presence 
of an emission line at E~6.88 keV in Obsl and an absorption 
line at E~7.75 keV in Obs2. We initially modelled these fea- 
tures with Gaussian lines. Their best-fitting parameters are 
reported in Table 2. The fit improvement with respect to the 
baseline models are Ax 2 =25.6/15.3 for 3/2 additional model 
parameters, for the emission/absorption line, respectively. 
Their detection significance is high, >99.99% and 99.95% 
from the F-test and >99.9% and 99.8% from 1000 Monte 
Carlo simulations, respectively. The Monte Carlo simula- 
tions were performed following the method of Tombesi et 
al. (2010a). We estimated the random detection probability 
for lines in the intervals E=6.5-7.5 keV and E=7-10 keV for 
the emission and absorption feature, respectively. 

We checked that the emission and absorption lines are 
independently detected in both XIS and XIS 3 and the 
results are consistent with the XIS 1. As reported in Ta- 
ble 1, the 4-10 keV background level is <1% of the source 
counts. However, we note that the XIS background has an 
instrumental emission line at E=7.5 keV due to Ni Ka (Ya- 
maguchi et al. 2006). This is not consistent with but close to 
the observed energy of the absorption line in Obs2 and we 
performed several tests to assure that the results are not af- 
fected by an erroneous background subtraction. We avoided 
regions too close to the chip corners, where the instrumen- 
tal contamination is higher, and checked that the results are 
consistent selecting different background regions. We also 
inspected that the observed energy of the background line 
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Table 1. Observation log and best-fitting baseline model parameters. 



Obs 


Date 


Exp 


Counts 


Flux 


r 


E 


I 


EW 


X> 


1 


2010/9/2 


59 


72/0.7 


2.97 ±0.02 


1.65 ±0.02 


6.38 ± 0.04 


2.3 ±0.3 


40 ±7 


1418.5/1380 


2 


2010/9/9 


59 


93/0.8 


3.78 ±0.02 


1.68 ±0.02 


6.39 ±0.02 


2.2 ±0.3 


33 ± 6 


1490.4/1469 


3 


2010/9/14 


52 


79/0.7 


3.68 ±0.02 


1.70 ±0.02 


6.40 ± 0.02 


2.7 ±0.4 


38 ± 10 


1427.0/1410 



Note. Columns: observation number; observation date; net XIS exposure in ks; 4-10 keV XIS 03 source/ background 
counts in units of 10 3 ; 4-10 keV flux in units of 10 -11 erg s _1 cm -2 ; photon index; Fe Ka rest-frame energy in 
units of keV; intensity in units of 10 -5 ph s _1 cm -2 ; equivalent width in units of eV; best-fit statistic. 



is not consistent at the 90% with that of the absorption 
line and the intensity is only ~30%. We checked that the 
absorption line parameters are consistent with and without 
background subtraction and that it is present also fitting the 
XIS and XfS 3 separately. Moreover, given that the back- 
ground level and instrumental emission line are identical for 
the three observations and also the source spectra are com- 
parable, if the absorption line was due to some systematics, 
it would have been present in all the spectra, and especially 
in Obsl which has the lowest flux. 

Only highly ionized Fe K emission lines are expected 
in the interval 6.5-7 keV, in particular Fe XXV Hea at 
6.7 keV and Fe XXVf Lya at 6.97 keV (Kallman et al. 2004). 
However, the energy of the Gaussian emission line in Obsl 
is E=6.83-6.93 keV at the 90% level and therefore is 
not consistent with these transitions, if not requiring an 
energy blue/red-shift. The emission line is resolved with 
cr=84 ± 38 eV. Therefore, we tested the possible modelling 
as a blend of rest-frame narrow Fe XXV-XXVI lines. How- 
ever, this provides a much worse fit (Ax 2 /A^=17/2) com- 
pared to that with a single Gaussian line profile. An alterna- 
tive and intriguing possibility is the identification with the 
blue peak of a relativistic line (e.g., Sambruna et al. 2011). 
Therefore, we replaced the emission line in Obsl with the 
relline model of Dauser et al. (2010). Given the limited 
S/N, the black hole spin, emissivity and outer radius can 
not be constrained and were fixed to the typical values of 
a=0, /3=— 3 and r o „t=1000 r g . We obtain comparable fits 
assuming in turn an energy of the line equal to the neutral 
Fe Ka, Fe XXV Hea or Fe XXVI Lya. The best fitting pa- 
rameters in the three cases are ri„< 13 r g and i=44° ± 2°, 
r in =9±2r g and i=32° ±2°, r in =57±l 3 6 r g and i=13°±4°, 
respectively. The fit improvement is Ax 2 /Az^=26/2, which 
corresponds to a confidence level >99.99%. However, if we 
require the disc inclination to be consistent with that of the 
jet of -18° (Jorstad et al. 2005) or 10° < i < 26° (Lewis 
et al. 2005), the solutions with ionized lines are preferred. 
The Fe XXVI model is shown in Fig. 2 (Panel 4). In partic- 
ular, if we adopt an inclination of i=18° and let the energy 
free to vary, we obtain E=6.94 ± 0.05 keV, which is con- 
sistent with Fe XXVI, and r in =20~tl° r g . The fit improve- 
ment is still high, Ax 2 /A^=25/2. The intensity of the line 
is 7=(2.4 ± 0.6) x 10~ 5 ph s" 1 cm" 2 and EW=54 ± 27 eV. 

The energy of the absorption line at 7.75 keV in Obs2 
is also not consistent with any known atomic transition. As- 
suming a conservative identification with the closest ones 
due to Fe XXV-XXVI, the line energy requires a blue-shifted 
velocity of ~0. 1-0. 14c. We modelled the absorber with the 
photo-ionization code Xstar (Kallman & Bautista 2001). 
We assumed a standard T=2 power-law continuum with 



cut-off at 100 keV and standard Solar abundances (Asplund 
et al. 2009). Due to the limited energy resolution of the 
XIS, the absorption line is not resolved, but we can place 
a 90% upper limit of 10,000 km/s. Therefore, we performed 
fits using three Xstar tables with turbulent velocities of 
1000 km/s, 3000 km/s and 5000 km/s and averaged the 
resultant fit parameters. To search for the possible best- 
fitting solution we stepped through the absorber redshift 
in small increments of Az=10~ 3 in the interval between 
0.1 to —0.4, leaving also the other parameters of the ab- 
sorber, i.e. Nh and log£, and the continuum free to vary. 
A slightly better fit is provided by the 3000 km/s table, 
which is shown in Fig. 2 (Panel 4). The best-fitting solu- 
tion yields iw=0.106±0.006c, log£=4.32±0.12 erg s _1 cm 
and N H = (7.7 ± 2.9) x 10 22 cm~ 2 . This is consistent with 
Fe XXVI being the dominant ionic species. The fit improve- 
ment is high, Ax 2 /A^=23/3, which corresponds to a detec- 
tion probability ^99.99%. 

The 90% upper/lower limits of the intensity and EW 
of the two additional emission and absorption lines in the 
three observations are reported in Table 2. A constant be- 
haviour of the lines can be tested by fixing the energy and 
width to the best-fitting values and comparing with the fits 
with free intensity/EW. We can rule out a constancy of the 
emission line in Obs2 and Obs3 at 99.7%. Instead, the con- 
stancy of the absorption line in Obsl and Obs3 is ruled out 
at 99.9% and 92%, respectively. Therefore, both features 
are significantly variable. Finally, consistent line parameters 
are found when the broad-band spectra are considered. The 
complete broad-band analysis and fits including ionized re- 
flection models, along with a comparison with previous X- 
ray (e.g., Tomesi et al. 2010a; Ballo et al. 2011) and radio 
observations of 3C 111 (e.g., Chatterjee et al. 2011) will be 
presented in forthcoming papers. 



3 DISCUSSION AND CONCLUSIONS 

From the spectral variability time-scales of ~7 days we 
can constrain the distance of the material(s) producing 
the emission/absorption features from the X-ray source to 
d < cAt < 1.8 x 10 16 cm (<0.006 pc), which corresponds 
to <40-600 r g assuming the black hole mass of M04 or 
Cll, respectively. This independently tells us that the ma- 
terial must be rather close to the SMBH. Phenomenologi- 
cally, the emission line at 6.88 keV in Obsl is in a region of 
the spectrum in which only emission lines from Fe XXV- 
XXVI would be expected, which corresponds to log£~3- 

4 erg s _1 cm (Kallman et al. 2004). The ionization parameter 
is defined as £,=L ion /nr 2 (Tarter et al. 1969), where Li on is 
the 1-1000 Ryd (1 Ryd=13.6 eV) absorption corrected ion- 



4 F. Tombesi et al. 



E I 





Obs1 


j 
















\ 


(4) 



Obs2 

. . .r . . . ( 1 >- 










(2) 




_v8v 






/ 


\ 








/ (4)_ 



; i 


' Obs3 

\ , : 

H , , i 1>- 








(» 










: / 






(4) 



Observed energy : KeVr 



Observed energy (key) 



Observed energy [fceV) 



Figure 2. Comparison of the XIS 03 spectra of the three observations zoomed in the 5-9 keV interval. Panel 1: ratios with respect to 
a Galactic absorbed power-law continuum. Data were binned to a S/N=30 only for plotting. Panel 2: A% 2 residuals. Panel 3: contour 
plots with respect to the baseline models including also the narrow Fe Ka emission line. The contours refer to Ax 2 levels of —2.3, —4.61 
and —9.21, which correspond to F-test confidence levels of 68% (red), 90% (green) and 99% (blue), respectively. The black indicates the 
reference level A\ 2 = +0.5. Panel 4 : final best-fitting models. Besides the baseline models they include a relativistic Fe XXVI emission 
line in Obsl and an Xstar photo-ionization grid with turbulent velocity of 3000 km/s in Obs2, respectively. See text for more details. 
The vertical lines indicate the observed energy of the narrow Fe Ka emission line. 



izing luminosity, n is the number density and r the distance 
from the source. Substituting L 4on ~5.8x 10 44 erg/s for Obsl 
extrapolated from the spectrum and the distance d, we ob- 
tain a density n>10 9 cm -3 . Given the compactness of the 
distribution and substituting in Nu^nd, we can estimate 
that the reflecting material is possibly Compton thick, with 
Nh>W 25 cm -2 . However, the measured line energy is not 
consistent with Fe XXV-XXVI at rest. An association with 
Fe XXV implies a blue-shift of ~0.026c. This could be ex- 
plained by the bulk of reflection from the approaching side 
of the disc or a large scale, Compton-thick outflow. However, 
the latter seems unlikely given that we do not observe a P- 
Cygni profile and simultaneous absorption lines. Instead, an 
identification with Fe XXVI implies a red-shift of ~0.013. 
If interpreted as a gravitational redshift, this indicates that 
the bulk of the emission is at <80 r g . Therefore, both phe- 
nomenological arguments and a detailed modelling with a 
relativistic profile suggest the line in Obsl is most probably 
generated by reflection off the accretion disc at ~20-100r fl 
from the SMBH. 

Concerning the blue-shifted absorption line in Obs2, 
the photo-ionization modelling indicates that the absorber 
is highly ionized, log£~4.3 erg s _1 cm, possibly mildly 
Compton-thick, Nh—8 x 10 22 cm -2 , and has a mildly rela- 
tivistic outflow velocity « otl t— 0.1c. This clearly suggests an 
association with a UFO (Tombesi et al. 2010a,b). Follow- 
ing similar reasonings to those in Tombesi et al. (2010a), 
we can estimate several physical parameters of the outflow. 
Throughout, we assume a covering fraction C~0.5 as esti- 
mated for the UFOs in Tombesi et al. (2010a, b). The ab- 
sorption corrected ionizing luminosity in Obs2 extrapolated 
from the spectrum is Li „~8.1 x 10 44 erg/s. The ~7 days 
variability indicates a compact absorber. Substituting the 
best-fit parameters in the £ definition, we obtain a distance 
from the ionizing source <5 x 10 17 cm, consistent with the 
variability upper limit d. On the other hand, substituting 
the distance d, we can estimate a density n>10 8 cm -3 . 



The bolometric luminosity can be derived from 
the relation Lt i — WLi 0n (McKernan et al. 2007). As- 
suming a constant flow, the mass outflow rate is 
M ut=4:TYC^ 2L mHVout r ^i Mq/ut and the mass accre- 
tion rate is M acc —Lhoi/Tic?^l M@/yr, using r/=0.1. Then, 
Mout/Macc ~1, suggesting that the instantaneous ejection 
of mass in Obs2 is comparable to that accreted by the SMBH 
in that moment. The ~40% increase in luminosity between 
Obsl and Obs2 could have been driven by an increase in ef- 
ficiency and/or accretion rate. Interestingly, from the ratio 
A/ °" f ",°" t ~1 we derive that the wind momentum is of the or- 
der of the radiation field or photon momentum and therefore 
radiation pressure could have played an important role in ac- 
celerating the wind (e.g., King & Pounds 2003; King 2010). 
The mechanical power of this outflow is EK = \Moutv 2 ut ^5x 
10 44 erg/s, which is comparable to the typical jet power 
for radio galaxies of ~10 43 -10 45 erg/s (Rawlings & Saun- 
ders 1991). Then EK^0.06Lt o i, again consistent with mo- 
mentum driven outflows (e.g., King & Pounds 2003; King 
2010). The Eddington luminosity and relative accretion rate 
are L Bdd ~2.6-39.0 x 10 46 erg/s and M E dd=L Edd /r 1 c 2 ^-~ 
50 Mq/ht, depending on the SMBH mass estimates of Cll 
or M04. The Eddington ratio L boi / 'L Ed d~M acc / M Ed d<0-3 
is moderate. However, the L Ed d should be regarded only as 
an upper limit for the luminosity needed to accelerate winds 
by radiation pressure because it takes into account only the 
opacity due to Thomson scattering. Considering only this 
term, we obtain a lower limit to the wind opacity in Obs2 
of t~ct t N h ~0.05. 

It is important to note that in Obs2 we observe the 
wind when it has already been accelerated at ~0.1c and, 
due to the low opacity and high ionization, it follows that 
it would be hard to further accelerate it through radiation 
pressure alone. However, the ratio of ~1 between the wind 
and photon momenta suggests that radiation pressure had 
an important role in launching the outflow. Given the high 
column density of ~10 23 cm -2 in Obs2, it is plausible that 
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during the launching and acceleration process the density 
was higher and therefore the wind was less ionized and pos- 
sibly optically thick. In this case, the many absorption lines 
and edges from light elements up to iron could have provided 
the additional opacity (e.g., Proga & Kallman 2004). Impos- 
ing r~l we can estimate the photospheric radius (King & 
Pounds 2003; Pounds et al. 2003), which indicates the region 
where the photon momentum was deposited into the wind, 
and we derive i? p h~100 r g . We note that the observed ve- 
locity of ~0.1c and column density are only conservative 
estimates, corresponding to the components along our line 
of sight. The first value corresponds to the escape velocity at 
^100-200 r g and, thus, the outflow is likely to leave the sys- 
tem even if does not undergo further acceleration. Moreover, 
the outflow velocity and ~7 days variability are comparable 
to the disc Keplerian velocity and dynamical time scale at 
that location. 

Therefore, we can derive a picture in which an increase 
in the X-ray illumination, possibly linked to a rise in effi- 
ciency/accretion rate, causes a truncation of the inner parts 
of the disc or overionizes the material responsible for the rel- 
ativistic line. This can explain the lack of the 6.88 keV emis- 
sion line in Obs2-3. Then, an outflow is lifted up from the 
disc surface at ~100r 9 and is accelerated by radiation pres- 
sure to a velocity of ~0.1c between Obsl-2. The outflowing 
plasma might then flow along the magnetic field lines and be 
further accelerated by magnetic torque and pressure (Kato 
et al. 2004; Oshuga et al. 2009; Fukumura et al. 2010). This 
latter consideration is plausible given the low inclination of 
the outflow with respect to the jet of ~ 18° and the fact that 
3C 111 is a well known superluminal source. In particular, 
this ionized outflow might be associated with external lay- 
ers of the jet and provide the required external pressure to 
collimate it (Kato et al. 2004; McKinney 2006). The possi- 
ble disappearence of the blue-shifted absorption line in Obs3 
could then be explained by an increase in ionization and/or 
the absorber moving away from the line of sight. The se- 
quence X-ray dip followed by ejection event is reminiscent 
of the disc-jet connection recently found in this source by 
Chatterjee et al. (2011). However, the actual link between 
the outflow and jet is not clear yet. If the corona, where 
the X-rays originate, is indeed related to the base of the jet, 
then the loading of electrons in the jet may cause a flux de- 
crease. Then the particles are released through an outflow 
and an increase in flux follows. In a subsequent paper we 
will report on a more detailed comparison with respect to 
the radio properties of this source and possibly the connec- 
tion between the accretion disc, outflows and the jet would 
be better quantified. 
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